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winch  irdkiqurt  will  punt  mint  uirtul  lor 
problem*  o\ihi*  size  fo  he  prepared  lm  the 
WOCE  hvd\graphK  data  vi.  11  is  proposed 
u>  nan  a  cooperative  utimuiutd  fllori  10  use 
i hr v  data  in  a\v*tematK  wav 
An  comparison  ol  different  existing  tech¬ 

niques  can  be  a  fVst  step  in  ihn  dim  lion  A 
number  ol  paitKijLms  weie  actively  interest¬ 

ed  in  hemg  imniiA  I  he  da  a  set  is  available 
in  alt  participants.  aVt  also  in  other  interest¬ 
ed  parties  (contact  (. \Wuntrh  11  'D  I.  »r  the 
mailing  hst  UtM'E  IxVeRSDE  on  telemat! 
f>MX‘F  T  (M'.FAN)  li  iVrxpcctnl  that  manv 
oilier  helds  will  hr  madras  attahle.  through 
<  cillaboraine  effort.  in  appropriate  Is  gruided 
lorni  these  should  incUnl\ transient  tracers. 

tar  lion  rel.nrci  helds.  imprtVrd  rsnmatrs  ol 

the  wind  held  during  the  pAiod.  and  surface 
infill  measurements  f  he  grrVip  agreed 
share  data  as  *»ell  as  model  relulis  and  code 
The  participation  ot  other  tcie^itu  is  wel¬ 
comed 
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Coupled  Ocean- 
Atmosphere  Models 

2  lit  International  Liege 
Colloquium  on  Ocean  Dynamics 

,'1/fin  I.  HVirufem 

Oflue  of  Njv  jl  Rr-vjnh  Ailinfion  Vj 

B.  Edward  M( Donald 

Naval  Oceanoeraph*  and  AtrnosphrrH  Research 
Laboratory.  Sirnmt  ^oace  Center.  Mai 

(>!i>bal  scale  mean  dsnamKi  related  10  El 
Niho- Southern  Oscillation  ttNSOi  was  the 
focus  of  the  2 1  si  International  Colloquium  on 
Ocean  DsnamKs  LXSO  provides  modelers  a 
great  deal  of  data  revealing  corretauon  be¬ 
tween  sea  surface  temperature  tSSTl  anoma¬ 
lies  and  atmospheric  pressure  differentials  in 
the  tropics  (See  tnf\rlA  1 1989|  for  an  excel- 
lem  review  of  the  data  I 

The  meeting  was  held  Mas  8-12.  1989.  ai 
the  University  ol  Liege.  Belgium  More  than 
50  scientists  attended,  from  Australia.  Bel 
gium.  Canada,  federal  Republic  ol  Carmans. 
Trance.  Hungary,  Italy .  |apan.  Netherlands. 
People  i  Republic  of  China,  Spain.  L.K..  L  S  . 

and  U  S  5  K. 

The  16th  Colloquium,  in  1984 

The  I9R9  colloquium  followed  one  on  the 
same  theme  »n  1984  To  give  contest  lor  the 
1989  meeting  it  is  helpful  to  summanre  the 
earlier  one  In  1981  the  ibth  International 
Liege  Colloquium  on  Ocean  Dvnamtcs  was 
the  first  meeting  to  be  desoted  to  coupled 
ocean-atmosphere  models  fSikoul.  19851  The 
motivation  lor  the  meeting  was.  m  large  part, 
the  vers  strong  F.I  Nino  of  1982-1985 

At  the  1984  colloquium.  43  scientists  pre¬ 
sented  papers.  Almost  SO  addressed  orran  re¬ 
sponse  to  speficted  atmospheric  forcing  or 
sice  versa.  Onlv  2  attempted  to  couple  ocean¬ 
ic  and  atmospherK  general  circulation  mocfels 
|C(Alt|  Thr  modeling  state-of  (he  art  at  that 
time  was  not  reads  to  attack  a  benchmark  on 
the  scale  of  the  ESSO  cscle.  One  problem 
♦hat  became  esident  in  coupled  W.Ms  was  a 
discouraging  downward  climate  drift  in  SST 
and  air  temperature.  On  the  whole,  the  held 
of  interactne  ocean-atmosphere  modeling  in 
1984  was  in  its  infancs.  despite  sears  of  expe¬ 
rience  with  separateis  developed  atmospheric 
and  oceanic  GCMs 

The  Slat  Colloquium,  in  1989 

At  the  1989  t  iege  colloquium  52  papers 
were  presented,  with  about  40*  addressing 
or  using  ocean -atmosphere  coupling  parame¬ 
ters.  R.  VV  Stewart  of  the  University  of  Vic¬ 
toria.  British  Columbia.  Canada,  opened  the 
meeting  on  the  optimistic  note  that  ocean 
models  are  now  approaching  the  point  where 
predictions  are  possible  He  said  that  political 
decision  makers  will  take  action  based  on  this 
predictive  capability  and  that  it  is  incumbent 
on  the  research  conn  lunns  to  express  the 
level  of  unceriamiv  that  exists  in  the  predic¬ 
tion*  As  the  conference  progressed,  it  he- 
came  clear  that  his  assessment  of  predictive 
capability  was  optimistic  and  his  warning  of 
uncertainty  well  founded 

The  papers  given  during  the  first  half  of 
the  meeting  dealt  primarily  with  responses  of 
one  medium  to  the  other,  with  little  concern 
for  coupled  interaction.  The  stated  theme  of 
ihe  meeting  was  coupling,  so  the  importance 
of  the  first  half  was  to  identify  the  relevant 
physics  and  numerics  rather  than  present 
specific  results. 

Scientist*  giving  papers  in  the  second  half 
of  the  meeting  reported  varsing  degrees  of 
success  in  modeling  the  ocean  and  atmo¬ 
sphere  together  .Some  findings  will  be  sited 
here  because  of  their  importance  to  the  cou¬ 
pling  issue  About  5  of  the  modeling  papers 
addressed  ENSO  directly  Some  promising 
sign*  in  the  model  results  appeared  in  the 
form  of  eastward^  propagating  SSI  anoma¬ 
lies  with  time  scales  ol  order  2  sears  lor 
crossing  the  Pacific  Ocean. 

Atmospheric  Dynamics 

Papers  on  the  first  das  dealt  largely  with 
atmospherK  response  to  El  XtAn.  Warm  SSI 
anomalies  have  two  immediate  effects  on  the 
atmosphere  above  them  increased  latent  heat 
Bus.  and  increased  convergence  Since  atmo¬ 
spheric  motion  is  in  sen  delicate  balance  in 
the  tropics  due  to  the  prevalence  of  ennvee- 
*ion  and  the  absence  of  gromnphk  control, 
ihe  warm  SST  anomalies  quickly  induce 
broad-scale  convergence  ol  surface  air  T  his 
convergence  intensifies  convective  updralts. 
causing  increased  condensation  and.  ultimate¬ 
ly.  precipitation 

Condensation  liberates  heat,  which  intensi¬ 
fies  the  surface  convergence.  Hence,  in¬ 
creased  precipitation  generally  follows  centers 
of  SST  warming  It  is  the  convergence  ol  sur¬ 
face  air  that  is  primarily  responsible  for  the 
increased  precipitation,  raiher  than  the  stm 
pie  mcr-sse  m  heat  flux  that  comes 

from  warmer  SST  t  Rismusson  ol  the  Uni¬ 


versity  of  Maryland.  College  Park  itrKr»t>rd 
available  data  on  f  NSO  and  stressed  the 
need  io  understand  phase  correlations  evi¬ 
dent  m  annual  biennial  and  lower  Irr.juem  * 
(4  to  5  vrar  periodsi  modes  of  atmospheric 
variability 

I  fie  mechanism  bv  which  equatorial  anno 
spherK  responses  propagate  from  the  iroptcs 
to  midlatitudes  is  not  clear  In  lact.  not  all 
model  Calculations  propagate  the  re*(>onse  at 
all  V\ hen  propagation  ocucs.  «  does  so  prt 
manlv  m  response  to  warm  anomalies  One 
suggested  propagation  mechanism  involves 
the  interaction  of  equaiortallv  induced  con¬ 
vergence  with  the  jet  stream  Ihe  suggestion 
is  that  the  equatorial  convergence  interrupts 
the  ronal  pattern  of  the  jei  stream,  increasing 
iti  wave  number  and  making  it  more  meridi¬ 
onal  in  naiure 

In  nme.  increased  equatorial  convection 
leads  to  increased  cloudiness  and  increased 
outgoing  radiation  Radiation  is  .he  least  un¬ 
derstood  and  most  complicated  heat-transler 
process  to  calculate  in  the  atmosphere,  so  u  is 
the  weakest  component  ol  atmospherK  mod¬ 
els  Since  radiation  plass  such  an  important 
role  in  tropKal  atmospheric  motion,  atmo¬ 
spheric  models  sulfer  in  thie  region  of  the 
world 

Fortunately,  the  other  atmospherK  process¬ 
es  have  been  simulated  numerically  lor  more 
than  1  decadei  in  numerical  weather  predic¬ 
tion  iNW'P)  models  for  ihe  midlaitiudes  Vir¬ 
tually  all  modern  atmospherK  models  used 
lor  large-scale  air  sea  interaction  research  use 
the  "pnmitise  equations  |PF>  of  atmospherK 
motion  Such  models  adequately  simulate  the 
relevant  physical  processes.  w»rh  the  excep¬ 
tion  of  radiation 

Ocean  Dynamic* 

In  contrail  to  a  long  history  w»th  atmo- 
iphenc  PF  models  ocean  model*  are  relative¬ 
ly  newer  I  he  l>4«  i  I  'xeao  model  currently 
in  use  for  research  was  developed  bv  8»vjh 
[  I9fi9|  A  maturity  of  the  conference  papers 
about  large-scale  ocean  simulation  used  one 
of  the  variants  ol  ihis  model  Other  models  tn 
use  make  physical  simpltcations  tvpihcd  bv 
approaches  such  as  geostrophv.  quaii-geos- 
trophv.  lumped  parameter  "box  ocean  ap¬ 
proaches.  and  one-or  rwo  laser  Laver  modeli 
relate  icxcan  evolution  to  horizontal  motion 
wuhin  vertically  homogenized  lavers.  with 
moving  lavrr  boundaries  describing  the  ocean 
surface  and  thermochne 

(Nean  models  respond  to  surface  winds 
and  to  evaporation  and  precipitation  differ¬ 
ence*  Wind  stress  provokes  the  quKkesi 
ocean  responses.  with  time  Kales  down  to  less 
than  an  hour— during  passage  of  a  storm,  for 
example  On  die  ;xher  extreme,  the  global 
thermohalme  circulation  operate*  on  a  much 
longer  time  Kale.  I  hi*  circulation  transports 
the  ocean  *  heat  hon/on  taffy  and  i*  driven  by 
density  anomalie*  that  can  be  due  to  differen¬ 
tial  heating,  evaporation  and  precipitation,  or 
vertical  mixing  The  glohai  thermohalme  cir¬ 
culation  time  is  on  a  decade  Kale 

Wind  increase  anomalies  have  two  effects 
on  the  ocean:  wind-driven  current*  and 
mixed-laver  deepening  Both  effect*  strongly 
influence  ocean  circulation  The  curl  of  wind 
stress  t»  the  principal  factor  tn  Urge-Kale 
wind  -driven  ocean  circulation  Mixed-laver 
drrjK-mng.  along  with  evapoiauon  and  pre¬ 
cipitation.  contributes  to  the  thermohalme 
circulation 

Fresh  water  input  tpretipnation  and  run- 
ofT)  and  evaporation  contribute  oppositely  to 
mixed-laver  processes  that  drive  large-scale 
circulation  cun  :mi  Evaporation  increases 
mixed-laver  salinii'  Because  saline  water  is 
more  dense  than  lre*h.  increased  salinity 
lead*  to  increased  convection,  mixed-laser 
deepening  and  SST  cooling  Conversely, 
whrn  fresh  water  input  exceeds  evaporation, 
the  ocean  freshens,  convection  ts  decreased 
and  the  mixed  laver  shallow*  Mixed  laser 
deepening  or  shallowing  set*  Up  a  horizontal 
pressure  gradient  that  manifests  tuelf  in  hori¬ 
zontal  currents,  this  process  is  the  principal 
component  of  the  thermohalme  circulation 

T wo  papers  described  ocean  models  used 
to  investigate  thermohalme  currents  tn  thr 
AtlaniK  and  the  Pacific,  respectively  Ihe  At 
laniK  ltuds  ptedKted  an  mlubtting  effect  on 
these  furtenl*  that  results  from  increased 
ouiput  of  the  St  Lawrence  River  Both  stud¬ 
ies  showed  that  the  existence  of  dual  steads 
states  was  plausible,  suggesting  the  possibility 
that  dillerent  emulation  patterns  mas  have 
existed  tn  the  past 

Coupled  Interaction* 

fdoAof  5ro/c  , Modeli  DtscUsuon  of  model 

coupling  brgan  on  the  third  day  of  the  meet 
mg  About  half  the  remaining  participants 
discussed  uue  coupling  Due  to  thr  compli¬ 
cated  nature  •»(  each  individual  mediums  re 
sponse  to  ihe  other,  truly  coupled  models  re 
quire  enormous  computet  resource*  or  con 
sidrraMe  simplification  ot  Kith  I  he  latter 
grrailt  compromises  reality  and  the  former  is 
very  expensive  Hence,  all  the  papers  on  true 
coupling  were  filled  with  introductory  qualih 
cations  such  as  cimplfic  anon  or  prrltmi 
naiy 

TOR  teen/  <m  p  I  <i 
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TOR  t rcml  ftfmi  p  1) 

In  (In*  quaMvmg  vein,  it  is  worthwhile 
highlighting  some  o(  these  preliminary  re¬ 
sults  relating  to  SSI  anomaly  propagation  S 
Philander  and  W  Lau  ol  the  Geophysical 
Fluid  DtnjmKi  Lalmratorv  tGFDL)  in 
Princeton.  N  |  .  found  that  the  direction  of 
SST  anomaly  propagation  wai  sensitive  10  the 
resolution  of  the  ocean  model  1-  Tokioka 
and  colleagues  at  the  Meteorological  Re¬ 
search  Institute,  in  I  sukuba.  Japan,  found 
eastward  propagating  SSt  anomalies  with 
coupled  5-le'el  atmosphere.  14  lf»el  ocean 
models,  both  having  a  common  horizontal 
grid  resolution  ol  V  *  I'  K  Brvan  ol  GFDL 
found  SSI  anomalies  slowlv  propagating 
eastward  in  the  Atlantic  Ocean  He  noted 
that  unless  a  proper  balance  of  moisture  flu* 
between  ocean  and  atmospheric  models  is 
maintained  the  Atlantic  s  heat  transport 
mechanism  is  impaired,  leading  to  downward 
climatological  drift  eventually  terminating  in 
an  ice  age 

It  appears  that  the  models  of  each  medium 
are  sufficiently  advanced,  the  interest  in  cou¬ 
pling  it  sufhoentlv  high,  and  available  com¬ 
puter  resources  ate  rapidly  enough  increas¬ 
ing  so  as  to  make  true  couplings  attainable 
But  the  present  preliminary  attempts  make 
the  findings  suit  tentative. 

Chmau  Onf'l  The  most  persistent  of  these 
earls  findings  was  that  both  models  exhibit 
climate  ^£i(i  toward  cooling  Ocean  and  at¬ 
mospheric  GCMs  using  the  PE  approach 
sohe  their  equations  (at  least  in  the  vertical 
direction)  bv  finite  difference  methods  Finite 
differences  impiioiiv  or  explicits  introduce 
numerical  diffusion  to  ensure  smoothness 
and  stabilits  of  the  solution.  Since  the  ocean 
has  a  large  tesersoir  of  cold  water  beneath  a 


re  (alive  Iv  shallow,  warm  surface  laser,  a  nu 
mciwaJ  procedure  that  m«es  this  deep  large 
rool  reservoir  wnh  the  shallow  warm  one 
tends  to  lower  the  surface  temperature  I  his 
artificial  diffusion  of  subsurface  water  into 
the  mixed  laser  is  a  serious  problem  lor  PF. 
ocean  models  The  atmospheric  radiation  bal¬ 
ance  problem  also  contributes  to  a  cold  chma- 
rologscal  drift  When  ocean  and  atmospheric 
models  are  coupled  this  drift  accelerates  and 
eventualls  weakens  the  coupled  tvstem  be- 
sond  realits 

Other  Issues 

Manx  topics  other  shits  those  mentioned 
here  were  discussed  at  the  1989  Liege  collo¬ 
quium  Two  papers  dealt  wnh  smallrr-scale 
model  coupling  one  on  quasi  steads  match¬ 
ing  of  ocean  and  atmospheric  boundary  las¬ 
ers.  and  one  on  hurricane  dynamics. 

Other  issues  toothed  on  during  the  confer¬ 
ence  were  polar  ice  cover,  global  atmospheric 
albedo,  and  the  greenhouse  effect.  Another 
potenualls  important  topic  not  covered  was 
the  time  and  space  Kile  mismatch  between 
atmospheric  and  ocean  processes.  During 
question  periods,  modelers  were  asked  for 
specifics  on  how  interface  conditMMis  are  han¬ 
dled  in  light  of  this  mismatch.  Answers  were 
short  comments  indicating  that  fluxes  are  in¬ 
terchanged  from  one  to  a  few  times  per  dav. 
Apparentls.  a  framework  for  consistent  flux 
exchange  has  vet  to  be  investigated.  It  seems 
likelv  that  the  temporal  and  spatial  Kale  mis¬ 
match  between  ocean  and  atmosphere  will  re¬ 
quire  special  numencaJ-lheoreticai  methods 
for  proper  treating  of  the  interface  The  rela¬ 
tive  importance  ol  these  issues  will,  one 
hopes,  be  clarified  bs  work  in  progress  dur¬ 
ing  the  next  few  vears 


Summary  and  Cone  I  ua  ions 

Similarities  and  differences  between  the 
Fust  and  second  Liege  colloquiums  descend  to 
ocean -atmosphere  coupling  show  the  pro¬ 
gress  that  has  been  made  in  the  hie  sears  fie 
tween  the  meetings  Although  both  were 
nominally  devoted  to  coupling,  both  had  a 
minority  of  papers  truls  devoted  to  that  top¬ 
ic  The  1989  colloquium  made  more  progress 
m  that  direction,  but  coupling  is  still  a  vers 
complicated  issue  and  most  investigators  prr 
fer  the  simpler  problem  of  one  medium  re¬ 
sponding  to  the  other  with  no  feedback 

Despite  this  preference,  a  number  of  indi 
cators  point  to  an  optimistic  direction  for 
coupled  modeling  In  198  4  the  community 
was  hardly  ready  to  attempt  senous  interac¬ 
tive  coupling  of  GCMs  because  of  limned 
computer  resources  and  few  proven  models 
Independent  ocean  and  atmospheric  model¬ 
ing  has  made  considerable  progress  over  the 
past  5  years  Computer  resources  are  consid¬ 
erably  greater  now  tfc;n  5  yea's  a*o.  and  ex¬ 
cellent  global  models  exist  for  coupled  ex¬ 
perimentation. 

A  few  anecdotes  help  capture  of  the  collo¬ 
quium.  Referring  to  the  history  ol  uncoupled 
models,  one  speaker  pointed  out  impropri¬ 
eties  ‘  of  "topless  oceans  and  bottomless  at¬ 
mospheres."  Another  speaker  compared  the 
ocean-atmosphere  interaction  to  a  marriage 
The  parlies  10  this  m  a  mage  exchange  heat 
and  have  periods  of  elevated  stress  At  times 
one  member  of  the  partnership  i«  slave  to  the 
other,  and  at  other  times  forces  the  other 
into  action  Each  member  of  the  partnership 
has  iu  own  personality:  the  atmosphere  is  im¬ 
petuous.  the  ocean  rather  reserved 

Another  speaker  pointed  cut  ;he  need  for 
reliable  component  models  with  a  profound 


observation  when  two  had  rrvxJels  are  <«»u 
pled.  H  is  possible  to  get  a  bail  resuli  \ 
more  serious  summary  was  given  bi  S  Inr 
son  of  the  Mciroroiogna]  Office  Bracknell. 

I  K  .  when  he  said  Coupled  modeling  u 
vers  difficult  but  we  are  making  progress 
Because  ocean  atmosphere  model  coupling 
is  a  relatively  new  and  uncharted  held  ambi¬ 
tious  modeling  efforts  are  soil  m  shake  down 
phases  Modelers  are  generally  moving  to¬ 
ward  qualitative  reaumabdiu  of  results  wnh 
quanuume  atcuraev  of  prediction*  set  to  he 
determined  In  comparison  wnh  develop¬ 
ments  in  the  component  modeling  areas  firm 
conclusions  and  new  insights  in  coupled 
ocean-atmosphere  modeling  are  elusive 
Wnhin  the  component  models,  phvsxs  im 
prosement  seems  to  be  required  in  atmo¬ 
spheric  radiation  simulations,  and  improved 
numerics  teems  to  be  required  to  remove  nu 
mencal  diffusion  from  the  ocean  modeli 
Larger  computers  will  be  required  to  allow 
both  models  to  be  run  in  sufficient  time  and 
space  resolution  to  include  all  the  relevant 
processes  Since  this  latter  suggesuon  contains 
an  element  of  brute  force,  thetr  is  an  tndt 
cation  of  need  for  more  theoretical  work  on 
consistent  model  coupling  allowing  for  rfitpa 
rate  time  and  space  vales  of  the  ocean  and 
atmosphere 
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